Using the Noether symmetry approach to probe the nature of dark energy 
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We propose to use a model-independent criterion based on first integrals of motion, due to Noether 
symmetries of the equations of motion, in order to classify the dark energy models in the context 
of scalar field (quintessence or phantom) FLRW cosmologies. In general, the Noether symmetries 
play an important role in physics because they can be used to simplify a given system of differential 
equations as well as to determine the integrability of the system. The Noether symmetries are 
computed for nine distinct accelerating cosmological scenarios that contain a homogeneous scalar 
field associated with different types of potentials. We verify that all the scalar field potentials, 
presented here, admit the trivial first integral namely energy conservation, as they should. We also 
find that the exponential potential inspired from scalar field cosmology, as well as some types of 
hyperbolic potentials, include extra Noether symmetries. This feature suggests that these potentials 
should be preferred along the hierarchy of scalar field potentials. Finally, using the latter potentials, 
in the framework of either quintessence or phantom scalar field cosmologies that contain also a 
non-relativistic matter (dark matter) component, we find that the main cosmological functions, 
such as the scale factor of the universe, the scalar field, the Hubble expansion rate and the metric 
of the FRLW space-time, are computed analytically. Interestingly, under specific circumstances 
the predictions of the exponential and hyperbolic scalar field models are equivalent to those of the 
ACDM model, as far as the global dynamics and the evolution of the scalar field are concerned. The 
present analysis suggests that our technique appears to be very competitive to other independent 
tests used to probe the functional form of a given potential and thus the associated nature of dark 
energy. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.+X 
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1. INTRODUCTION 

The analysis of the available high quality cosmological 
data (supernovae type la, CMB, galaxy clustering, power 
spectrum etc) have converged during the last decade to- 
wards a cosmic expansion history that involves a spatially 
flat geometry and a recently initiated accelerated expan- 
sion of the universe (see [l|-|g| and references therein). 
From a theoretical point of view, an easy way to explain 
this expansion is to consider an additional energy com- 
ponent, usually called dark energy (DE) with negative 
pressure, that dominates the universe at late times. The 
simplest DE candidate corresponds to the cosmological 
constant (see 0-11 for reviews). Indeed the so called spa- 
tially flat concordance ACDM model, which includes cold 
dark matter (DM) and a cosmological constant, (A) fits 
accurately the current observational data and thus it is 
an excellent candidate to be the model which describes 
the observed universe. 

However, the concordance model suffers, among other 
[Icj . from two fundamental problems: (a) The "old" 
cosmological constant problem (or fine tuning problem) 
i.e., the fact that the observed value of the vacuum en- 
ergy density (p\ — Ac 2 /SttG ~ ICU 47 GeV A ) is many 
orders of magnitude below the value found using quan- 
tum field theory Q, and (b) the coincidence problem [TTJ 
i.e., the fact that the matter energy density and the vac- 
uum energy density are of the same order (just prior to 
the present epoch), despite the fact that the former is 



a rapidly decreasing function of time while the latter is 
stationary. 

Attempts to solve the coincidence pro blem have been 
presented in the literature (see [1, [9l Il2|| and references 
therein), in which an easy way to overpass the coinci- 
dence problem is to replace the constant vacuum en- 
ergy with a DE that evolves with time. Nowadays, the 
physics of DE is considered one of the most fundamental 
and challenging problems on the interface uniting As- 
tronomy, Cosmology and Particle Physics. In the last 
decade there have been theoretical debates among the 
cosmologists regarding the nature of this exotic compo- 
nent. Many candidates have been proposed in the lit- 
erature, such as a cosmological constant A (vacuum), 
time-varying A(t) cosmologies, quintessence, k— essence, 
quartessence, vector fields, phantom, gravitational mat- 
ter creation, tachyons, modifications of gravity, Chaply- 
gin gas and the list goes on (see @, fl3H29j and references 
therein). 

In the original scalar field models [3l| and later in the 
quintessence context, one can ad- hoc introduce an adjust- 
ing or tracker scalar field <j> [3 (different from the usual 
SM Higgs field), rolling down the potential energy V((f>), 
which could resemble the DE [8|, |{§ I32l435j . However, it 
was realized that the idea of a scalar field rolling down 
some suitable potential does not really solve the problem 
because 4> has to be some high energy field of a Grand 
Unified Theory (GUT), and this leads to an unnaturally 
small value of its mass, which is beyond all conceivable 



2 



standards in Particle Physics. As an example, utilizing 
the simplest form for the potential of the scalar field, 
V{4>) = rr& <fi 2 /2, the present value of the associated vac- 
uum energy density is p\ — (V ((/>)) ~ 10~ n eV 4 , so for 
(</>) of order of a typical GUT scale near the Planck mass, 
Mp ~ 10 19 GeV, the corresponding mass of 4> is expected 
in the ballpark of ~ H ~ 1CT 33 eV. 

Notice that the presence of such a tiny mass scale in 
scalar field models of DE is generally expected also on the 
basis of structure formation arguments (36l - f38| ; namely 
from the fact that the DE perturbations seem to play 
an insignificant role in structure formation for scales well 
below the sound horizon. The main reason for this homo- 
geneity of the DE is the flatness of the potential, which 
is necessary to produce a cosmic acceleration. Being the 
mass associated to the scalar field fluctuation propor- 
tional to the second derivative of the potential itself, it 
follows that 7710 will be very small and one expects that 
the magnitude of DE fluctuations induced by (j> should be 
appreciable only on length scales of the order of the hori- 
zon. Thus, equating the spatial scale of these fluctuations 
to the Compton wavelength of <p (hence to the inverse of 
its mass) it follows once more that < Hq ~ 10~ 33 eV. 

Despite the above difficulties there is a class of viable 
models of quintessence based on supersymetry, super- 
gravity and string-theory which can protect, under of 
specific potentials, the light mass of quintessence (for a 
review see 39] and references therein). In spite of that, 
this class of DE models have been widely used in the lit- 
erature due to their simplicity. Notice that DE models 
with a canonical kinetic term have a dark energy EoS 
parameter —1 < w<f, < —1/3. Models with (w^ < —1), 
sometimes called phantom DE [jo| . are endowed with a 
very exotic nature, like a scalar field with negative ki- 
netic energy. In any case, in order to investigate the 
overall dynamics of the universe we need to define the 
functional form of the potential energy. As we have al- 
ready mentioned the issue of the potential energy has 
a long history in scalar field cosmology and indeed sev- 
eral parametrizations have been proposed (exponential, 
power law, hyperbolic etc). 

The aim of the present work is to investigate which 
of the available scalar field potentials can accommo- 
date basic geometrical symmetries (connected to the 
space-time) namely Lie point and Noether. In fact 
the idea to use Noether symmetries as a cosmologi- 
cal tool is not new. In particular, it has been pro- 
posed that the existence of such symmetries are re- 
lated with conserved quantities and thus they can be 
used as a selection criterion in order to discriminate the 
dark energy models, includin g th ose of f(R) gravity (see 
[I3,0,H,0>H>|Ii>|43>H)- From a mathemati- 
cal point of view, the Lie point/Noether symmetries play 
a vital role in physical problems because they provide 
Noether (first) integrals, which can be used in order to 
simplify a given system of differential equations and to 
determine the integrability of the system. A fundamen- 
tal approach to derive the Lie point and Noether sym- 



metries for a given dynamical problem moving in a Rie- 
mannian space has been proposed recently by Tsamparlis 
& Paliathanasis I49B (a similar analysis can be found in 

The structure of the article is as follows. The basic 
theoretical elements of the problem are presented in sec- 
tion 2, where we also introduce the basic FLRW cos- 
mological equations for various potentials of the scalar 
field. The geometrical symmetries of the scalar fields and 
their connections to the potential energy are discussed in 
section 3. In section 4 we provide for a first time (to 
our knowledge) analytical solutions in the light of either 
quintessence or phantom scalar field cosmologies that in- 
clude non-relativistic matter (dark matter). Finally, the 
main conclusions are summarized in section 5. 



2. COSMOLOGY WITH A SCALAR FIELD 

The scalar field contribution to the curvature of space- 
time can be absorbed in Einstein's field equations as fol- 
lows: 
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k = 8ttG 



(1) 



where R^ v is the Ricci tensor and T Ml/ is the total energy 
momentum tensor T^ v given by = + T^cj)). 
Here T^ u (<f>) is the energy-momentum tensor associated 
with the scalar field <fi, and T M „ is the ordinary energy- 
momentum tensor of matter and radiation. Modeling 
the expanding universe as a perfect fluid that includes 
radiation, matter and DE with 4— velocity Uu, we have 
T MV = -Pg^u + {fl + PW^Uv, where p = p m + p<j, and 
P = Pm + P<p are the total energy density and pressure of 
the cosmic fluid respectively. Note that p m is the proper 
isotropic density of matter-radiation, p^ denotes the den- 
sity of the scalar field and P m , P$ are the corresponding 
pressures. In the context of a FLRW metric with Carte- 
sian coordinates 

ds 2 = -dt 2 + a 2 (t)- ^—~ 2 {dx 2 + dy 2 + dz 2 ) (2) 

(1 + 4' X ) 



the Einstein's field equations ([lj , for comoving observers 
(If 1 = <5q), provide 
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(3) 
(4) 



where the over-dot denotes derivative with respect to the 
cosmic time t, a(t) is the scale factor of the universe, and 
K-3 = 0, ±1 is the spatial curvature parameter. Also, the 
contraction of the Ricci tensor provides the Ricci scalar 



R = gTRau = 6 



(5) 



3 



Finally, the gravitational field equations boil down to 
Friedmann's equation 



H 2 = 



k 

= 3 (Pm + P4>) 



K3 

a 2 : 



and 



3H 2 + 2H 



-k( P rr 



P*)-^- 



(6) 



(7) 



where H{t) = a/a is the Hubble function. The Bianchi 
identity (which insure the covariance of the theory) 
T^v — amounts to the following generalized local 
conservation law: 



field depends only on the cosmic time t, and consequently 
4> tV = 0(5° where 4> — 4?. A scalar field 4>(t) with a po- 
tential V(4>) is defined by the energy momentum tensor 
of the form (for review see [39| and references therein) 



T^(<f>) 



2 S(^L^) 



(14) 



where is the Lagrangian of the scalar field. Although 
in the current analysis we study generically, as much as 
possible, the problem we will focus on a scalar field with 



L 4> = - V -^ V ^^, V -V{4>) 



(15) 



or 



p m + P4> + 3H{ Pm + P m + P<j> + P^) = . (8) 

Note that the latter quantities obey the following rela- 
tions: 



( Pm ,P m ) = (-Iff, 2?) ( P4> ,P^ = (-Iff (0), 3j(0)) 



Combining cqs. ([5]), © and ((HJ we obtain 
k , 
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(10) 



Assuming negligible interaction between matter and 
scalar field, eq. <j8j) leads to the following independent 
differential equations 



p m + 3H(p m + P rn ) = , 

P<j> + 3H(p^ + P ) = . 



(11) 
(12) 



In this work, we will present the global dynamics of 
the universe in the presence of a barotropic cosmic fluid 
whose the corresponding EoS parameters are given by 
w m = Pm/ Pm an d Wtf, — P^/p^. In what follows 
we assume a constant w m which implies that p m — 



PmO a 



-3(l+«)„ 



(cold w m = and relativistic w m = 1/3 
matter), where p m o is the matter density at the present 
time. Generically, some high energy field theories sug- 
gest that the dark energy EoS par ameter is a function of 
cosmic time (see, for instance, [30]) and thus 



P<f,o exp 



1 3[1 + WlP {a)} 



(1(7 



(13) 



where p^o is the DE density at the current epoch. 



2.1. The scalar field 

We consider a scalar field in the FRLW cosmology 
which is minimally coupled to gravity, such that the field 
satisfies the Cosmological Principle that is, <f> inherits the 
symmetries of the metric. This means that the scalar 



where 



H = o4 2 - v{4>) 



1 Quintessence 
— 1 Phantom . 



(16) 



(17) 



Therefore, using the second equality of eq.©, cq. (fl"4)) 
and eq. (|T6|) the energy density p^ and the pressure P$ of 
the scalar field are given by 
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and 



P^Tl^)=L^-etf-V{cp) 



(18) 



(19) 



Inserting eq. (|T8)) and cq. (|T9|) into eq.(fl2j) it is routine 
to derive the Klein-Gordon equation which describes the 
time evolution of the scalar field. This is 



> + -acj) + eVrf = 



(20) 



where — dV/dcf>. Obviously, if we use the current 
functional form of then eq.([7]) takes the form: 



a 
a 
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(21) 



The corresponding dark energy EoS parameter (de- 
fined before) is 



W0 



v(4>) 



P4> e(cf> 2 /2) + V{<f>) 



(22) 



The quintessence (e = 1) cosmological model accom- 
modates a late time cosmic acceleration in the case of 
Wtj, < —1/3 which implies that cp 2 < V((f>). On the other 
hand, if the kinetic term of the scalar field is negligible 
with respect to the potential energy [£- <C V{4>)] then 
the equation of state parameter is ~ — 1. In the case 
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of a phantom DE (e = —1), due to the negative kinetic 
term, one has < — 1 for (<ft 2 /2) < V{<p). 

The unknown quantities of the problem are a(t), (f>(t) 
and V(4>) but we have only two independent differential 
equations available namely eqs. (|20|) and (|2Tj). Thus, 
in order to solve this system of differential equations we 
need to assume a functional form of the scalar held poten- 
tial energy, V(cj>). In the literature, due to the unknown 
nature of the DE, there are many forms of potentials (for 
a review see (39|) which describe differently the physical 
features of the scalar held. Let us now briefly present var- 
ious potentials whose free parameters can be constrained 
by using the current cosmological data. 



• The power law potential [13j, [lj 

M i+n 



• The exponential potential [58j: 

V{<j>) = V exp(-Vk \cf>) 



(23) 



(24) 



• The unified dark matter potential (hereafter UDM) 
591: 



V{<j>) = V [1 + cosh 2 (A0)] 



(25) 



• The Pseudo-Nambu Goldstone Boson potential 

V{<f>) =//[! + cos(0//)] (26) 

• The exponential with inverse power [6lj |: 

V{4>) = M[exp( 7 /0) - 1] (27) 

• The supergravity motivated potential jl9j |: 

V{4>)= Mexp(0 2 )/</> 7 (28) 

• The early dark energy potential [62| : 

V{4>) = Vi exp(-Vfc Xi4>) + V 2 exp(-Vfc \ 2 4>) . (29) 

Interestingly the potential: 

V(</>) = Vb[cosh(A0) - l] p 

provides predictions which are similar to those of 
the early DE model, as far as the global dynamics 
of the universe is concerned [631 . 



• The Albrecht-Skordis model [64f : 

Vf» = M[A + {(j> - Bf] exp(- 7 0) 



(30) 



• The Chaplygin gas from the ordinary scalar field 
viewpoint |39l |: 



(31) 



Detailed analysis of these potentials exist in the litera- 
ture, including their confrontation with the data (see [39| 
for extensive reviews). It is worth pointing out that for 
some special cases attempts to find analytical solutions 
can be found in [6l|, l65T - [7l| (and references therein) . 



3. DARK ENERGY VERSUS SPACE-TIME 
SYMMETRIES 

In the last decade, a large number of experiments have 
been proposed in order to constrain DE and study its 
evolution. Naturally, in order to establish the evolu- 
tion of the DE equation of state a realistic form of H(a) 
is required while the included free parameters must be 
constrained through a combination of independent DE 
probes (for example SNIa, BAOs, CMB etc). However, 
there is always a range of parameters of the considered 
scalar held cosmological models for which a good fit with 
the observational data is provided. This implies that 
such DE models can not be distinguished observationally, 
since they provide similarly evolving Hubble functions. 

Practically, the goal here is to dchnc a method (selec- 
tion criterion) that can distinguish the DE models on a 
more fundamental (eg. geometrical) level. According to 
the theory of general relativity, the space-time symme- 
tries (Killing and homothetic vectors) via the Einstein's 
field equations (see eqfT]) , are also symmetries of the en- 
ergy momentum tensor (the matter generates the gravi- 
tational field). Owing to the fact that the scalar held is 
minimally coupled to gravity and it evolves in space-time 
one would expect that the scalar held must inherit the 
symmetries of the space-time as gravity does. 

It is interesting to mention that besides the geometric 
symmetries one has to consider the dynamical symme- 
tries, which are the symmetries of the field equations. 
These latter symmetries are known as Lie symmetries. 
In case the held equations are derived from a Lagrangian 
there is a special class of Lie symmetries, the Noether 
symmetries, which have the characteristic that lead to 
conserved currents or, equivalently, to first integrals of 
the equations of motion. The Noether integrals are used 
to reduce the order of the field equations or even to solve 
them. Therefore a sound requirement, which is possible 
to be made in Lagrangian theories is that they admit 
extra Noether symmetries. This assumption is model 
independent, because it is imposed after the held equa- 
tions have been derived, therefore it does not lead to con- 
flict whereas with the geometric symmetries while, at the 
same time, serves the original purpose of a selection rule. 
Since the basic equations in the scalar field cosmologies 
follow from a Lagrangian we can apply the above ideas 
by looking for scalar field cosmologies which admit extra 
Noether symmetries. 

In particular, let us consider a cosmological model 
which accommodates a late time accelerated expansion 
and it contains a scalar field <fr(t), described by a po- 
tential V((f>). We pose the following question: For the 
scalar field that lives into a 2- dimensional Riemannian 
space {a, (f>} and which is embedded in the space-time, 
how many (if any) of the previously presented potentials 
(see section 2.1) can provide non trivial Noether sym- 
metries (or first integrals of motion)! As an example, 
if we find a cosmological model (or a family of models) 
for which its scalar held produces non-trivial number of 
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first integrals of motion with respect to the other DE 
cosmological models, then obviously this model contains 
an extra geometrical feature. Thus, we can use this ge- 
ometrical characteristic in order to classify the explored 
DE cosmolo gica l model into a distinct category (see also 
mj3,H3,Iil)- Below, we present the geo- 
metrical method used in order to define the basic sym- 
metries, namely Lie point and Noether which lead to first 
integrals of motion. 

3.1. Lie and Noether symmetries 

We briefly present the main points of the method used 
to classify all two dimensional Newtonian dynamical sys- 
tems, which admit Lie point/Noether symmetries 1 . A 
first important ingredient is the use of two theorems 
which relate the Lie point and the Noether symmetries of 
a dynamical system moving in a Riemannian space with 
the special projective group and the homothetic group 
generators of the space respectively. These theorems have 
been given by some of us in [49} . In particular in a recent 
paper Tsamparlis & Paliathanasis [49| have provided an 
alternative way to solve the system of Lie point/Noether 
symmetry conditions (see their tables 1-15), for second 
order equations of the form: 

x l + T) k x j x k =F\ (32) 

Here T^ k (x r ) are general functions, along the solution 
curves and F % {x^) is a C°° vector field. Basically, equa- 
tions (|32j) are the equations of motion of a dynamical 
system in a Riemannian space in which the functions 
Tj k (x r ) are the connection coefficients (Christofell sym- 
bols) of the metric of the space (in our case {a, <p} see 
below). The key point, (see [49|), is to express the sys- 
tem of Lie point/Noether symmetry conditions of ea. (|32[) 
in terms of collineation (usually referred as symmetries) 
conditions of the metric. If this is achieved, then the Lie 
point/Noether symmetries of eo. (|3"l2)) will be related to 
the collineations of the metric. Therefore the determina- 
tion of the Lie point/Noether symmetries of eq. (|32p will 
be transferred to the geometric problem of determining 
the generators of a specific type of collineations of the 
metric. Then it will be possible to use the plethora of re- 
sults of Differential Geometry on collineations to produce 
the solution of the Lie point/Noether symmetry problem. 

The natural question to ask is: How one will select the 
Lie point/Noether symmetries of two different dynami- 
cal systems, which move in the same Riemannian space? 
The answer is simple. The left hand side of eg . (|3"!2)) con- 
tains the metric and its derivatives and it is common to 



Note that the Noether symmetries are a sub-algebra of the alge- 
bra defined by the Lie symmetries [49ll . A dot over a symbol in 
eq. (1-S21 indicates derivation with respect to the parameter s (the 
affine parameter along the trajectory - cosmic time in our case) 



all dynamical systems moving in the same Riemannian 
space. Therefore geometry enters in the left hand side 
of eq. (|32|) only. Each dynamical system is defined by 
the force field F l , which enters into the right hand side 
of ea. (|32l) only. Therefore, there must exist constraints, 
which involve the components of the Lie point/Noether 
symmetry vectors and the force field, and which will have 
to be satisfied for a collineation of the metric in order for 
it to be a Lie point/Noether symmetry of the specific 
dynamical system. A similar approach can be found in 

[El,H0&@,[H r -i 

In a subsequent work [49( Tsamparlis & Paliathanasis 
determined (among others) all two dimensional poten- 
tials which admit at least one Lie point and/or Noether 
symmetry. The results of the calculations have been col- 
lected for convenience in tables 1-15. From these tables 
one can read directly the aforementioned potentials and, 
furthermore, for each potential the admitted Lie point 
and Noether symmetries together with the correspond- 
ing Noether integrals. It is emphasized that no extra 
calculations are required. In the following we shall make 
use of these results. 

3.2. Using the geometry of the space {a, <f>} to 
constrain Dark Energy 

In this section we apply the previous theory into the 
scalar field cosmology. Interestingly, we can easily prove 
that the main field equations (l20l) and (|2"Tj) . described 
in section 2, can be produced by the following general 
Lagrangian: 

L = -3aa 2 + fca 3 L + ka 3 P m + 3K 3 a (33) 

in the space of the variables {a, (f)}. Also the action is 

S = J Ld 3 xdt . (34) 

Therefore, utilizing eq. dTtjl) and eq.(p?3"|) we can also obtain 
the Hamiltonian of this system 

E = \ {~ 6ah2 + ka ^ e ^) + ka3 1^0) - p m] - 3K 3 a . 

(35) 

In order to compute the Lie point/Noether symmetries 
of equations of motion ([20]) and (|2"Tj) . we consider the 
Lagrangian as the sum of a kinetic energy which defines 
the metric in the space of {a, </>} and an external force 
field. In particular, this 2 dimensional metric takes the 
form 

ds 2 = -6ada 2 + ka 3 ed<p 2 (36) 
which implies that 

1 <f> 3 k 

Obviously, in the case of phantom cosmology e = — 1, 
eq. (1361) points that that we have to replace <f> with i(f> for 
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mathematical convenience. Using now eq. (|37p and in- 
serting the variables x l — a(t), <j>(t) variables into eq. (j52")) 
we find that 



U+ Ya a 



1 -2 , * 



pa 



, + -d<£ = F*. 



Comparing with the equations of motion PU|) and ([2"T]) . 
we can define the external "forces" in terms of the scalar 
field potential V((j>) 



F a = ~K 3 -^[P m -V(<t>)} , F* 



-eV, 



On the other hand using the above Tj k functions we 
find after some simple algebra that the curvature of the 
{a, (/)} space is R — implying flatness (all 2 dimen- 
sional spaces are Einstein spaces hence R = implies 
the space is flat). Also, the signature of the metric 
eq. (l37T) is —1, hence the space is the 2-d Minkowski space. 
Therefore, according to theorems 1 and 2 of [49| the 
Lie point/Noether symmetries of the equations (I2TH) and 
(|2~T|) follow from the special projective group of the 2-d 
Minkowski space. To find explicitly these vectors and 
thus the corresponding Noether symmetries (or first in- 
tegrals), we have to bring the 2-d metric of eg. (1551) into 
its canonical form (i.e. ds 2 = —dx 2 + dy 2 ). Changing 
now the variables from (a, </>) to (r, 9) via the relations: 



the 2 dimensional metric ([M]) is given by 
ds 2 = -dr 2 + r 2 d0 2 



(38) 



(39) 



that is, (r, 9) are hyperbolic spherical coordinates in the 2 
dimensional Minkowski space {a, <p\. Next we introduce 
the new coordinates {x, y) with the transformation: 



x — r cosh (0) 
y = r sinh (9) 

which implies that eq. (|3"9")) becomes ds 2 = —dx 2 
We also point here that 

r 2 = x 2 — y 2 9 — arctanh(y/a;) . 

The scale factor (a(t) > 0) is now given by: 



3(z 2 - y 2 ) 



1/3 



(40) 
dy 2 . 

(41) 
(42) 



which means that the new variables have to satisfy the 
following inequality: x > \y\. 

In the new coordinate system (x, y) the Lagrangian 
and the Hamiltonian (1351) are written: 



(if -x 2 )- V e ff(x,y) 



(43) 



E=±(y 2 -x 2 )+V ef f(x,y) 



(44) 



where 

V eff (x,y) = (x 2 -y 2 ) 
Note that we have used 

i 

3 



Pn 



SK3 [x 2 — y 2 ) 3 . 

(45) 



3fc 



K a =Kal-\ Pm = —P,r 



and 



V(0) = ™V(9) 



(46) 



(47) 



If the matter pressure P m is constant then the Lagrangian 
(or Hamiltonian) is time independent, thus the system is 
autonomous implying that theorems 1 and 2 of [49j apply. 
In this case we also have (for more details see appendix 
A) that p m = j^g- — P m which obeys eq. (ITU . 

We now proceed in an attempt to provide the Lie point 
and Noether symmetries of the current dynamical prob- 
lem. Note, that for simplicity in the analytical treatment 
below, unless explicitly stated, we consider spatially flat 
FRLW (K3 = 0) scalar field models that include non- 
relativistic matter (P rn = with p m = p). Taking 
the latter into account the effective potential of eq. (|4"5")) 
is given by 



V e ff(x,y) = {x 2 



X 



r 2 V(9) 



(48) 



For this effective potential and with the aid of [5CJ and 
[49l | we find the following cases that admit Lie point and 
Noether symmetries: 

• Trivial Noether symmetries: First of all, it 
is well known (see for example [49| and refer- 
ences therein) that for a general effective potential 
V e ff{x,y) we have only the Lie symmetry dt- Ad- 
ditionally, for effective potentials given by eq.(| 
we have that: 

V eff (x,y)=x 2 (l-t I )v(l 
x z \x 

hence the case of table 8, line 5 (with d — 0) of [49[ 
applies and we have the additional Lie symmetry 
xd x + yd y . A linear combination between dt and 
xd x + ydy provides the following Lie symmetry: 



Xr 



c\d t + c 2 (xd x + ydy) 



(49) 



From the corresponding tables of theorem 2 of [49| , 
we see that in this case there is only the trivial 
Noether symmetry 2 d t , whose Noether integral is 



2 The operator d z denotes d/dz. Also, there is a difference of 
the results presented in |49( due to the Lorentzian character of 
the metric. This affects only the rotational part of the metric 
(non gradient Killing vectors) and gives ydx + xd y instead of the 
Euclidean yd x — xd y . 
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the Hamiltonian E =constant (dtE — 0). As it is 
expected, this result implies that all the scalar field 
potentials described in section 2.1 admit the trivial 
Nocther symmetry, namely energy conservation, as 
they should. 

• Extra Noether symmetries: Now we are look- 
ing for first integrals beyond the standard one. 
In particular, we are interested to check whether 
the scalar field potentials mentioned in this pa- 
per (see section 2.1) can admit non-trivial Lie 
point/Noether symmetries. The agrument is sim- 
ple: if for a given potential we find extra Noether 
symmetries which are related with conserved quan- 
tities then this particular model has an enhanced 
physical meaning (see also [Hj|, , [H , [11] ) • From 
the mathematical point of view the existence of ex- 
tra integrals of motion points the existence of an 
analytical solution (see next section). The novelty 
in the current work is that we find that the expo- 
nential potential (see ea IM|) and the UDM potential 
(see ea l2"5")) can be clearly distinguished from the 
other dark energy potentials because these are the 
only potentials from the list presented in section 
2.1, that accommodate extra Noether symmetries. 

Hyperbolic - UDM Potential: Generically, we 
use the following potential: 



V{8) 



coi cosh 2 (8) — L02 sinh 2 (9) 



or 



V eff (x,y) = r 2 V(9) 



UJlX 



u 2 y 



(50) 



(51) 



where we have used eqs. (|40l) and (|481) . The corre- 
sponding Noether symmetries, X n , are known (see 
for example [50j). These are (for wiw 2 ^ other- 
wise see appendix B): 

X ni = d t , X m = sinh (y/uJ^t) d x , X„ 3 = cosh (y/unt) d x 
X ni = sinh (-y/uit) d y , X ns = cosh (y/u^t) d y 

The Noether integrals are the Hamiltonian and the 
quantities: 



I n2 = sinh (y/Uit) x — 

l ll3 = COsh (y/uJlt) X - 

I nA = sinh (y/UJ^t) y - 
7„ 5 = cosh (y/uJ^t) y - 



Obviously the UDM potential is a particular case 
of the current general hyperbolic potential. Indeed 
for uji = 2uj 2 and with the aid of eas. ([55)l . we 
fully recover the UDM potential (see section 2.1) 



wi cosh (y/Uit) x 
Ui sinh (y/Uit) x 
uJ^ cosh (y/uJ^t) y 
cJJsinh (y/U^t) y 



where Vq = modulus a constant. 

Exponential Potential: Here we provide for a 
first time (to our knowledge) the Lie point and the 
Noether symmetries of the exponential potential. 
Indeed from table 8, line 3 (with d ^ 0) of [4!| one 
can immediately see that 

V eff (r,6)=r 2 V(e)=r 2 e- de . 

The corresponding Lie symmetry vector is: 



X L = 2td t + - {yd x 
a 



xd y ) 



(53) 



Concerning the Noether symmetries from table 14 
of we find that the Noether symmetry of the 
system for the potential V{9) = e~ dS is: 



X n = 2td t 



7> 1 *' ■' 



y- 



-x I i)„ 



(54) 



In general the Noether inte gral for the vector X n 
2td t + rfdi is (see eq.66 of [49j|): 



I = 2tE - rfgiji 1 

where -q 1 = (x + |y) d x + (y + d v 
algebra we compute: 



(55) 
After some 



'I'fhj'l'y = {x+^y y + \x ) 



-1 
1 



i.e. 



= - { x + -y J x + [y + -x ) y. 



21 E [x + -yjx- (y+2 x )fi 



(56) 



where E is the Hamiltonian. Using V — e~ to- 
gether with eq. (|3"8")) and ea P?)) we write the poten- 
tial to its nominal form (see section 2.1). This is 



V(<f>) = U exp 




(57) 



where Vq 



^ modulus a constant 3 . 



From now on, we focus on the exponential and the 
UDM potentials because they contain non trivial inte- 
grals of motion, implying the existence of exact analytical 
solutions (see next section). In section 4 we provide for 
a first time (to our knowledge) such analytical solutions 
in the light of either quintessence or phantom scalar field 
cosmologies that include also a non-relativistic matter 
(cold dark matter) component. 



v{4>) = v Q 



1 + cosh z — 



(52) 



3 In the special case of d = 2, the system admits an additional Lie 
symmetry d x + d y , with Noether integral I = x — y. 
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4. ANALYTICAL SOLUTIONS IN THE FLAT 
SCALAR FIELD COSMOLOGY 

In this section, we proceed in an attempt to analyti- 
cally solve the differential eqs. (|2"0)) and (f2"Tj) . We remind 
the reader that this is possible because the dynamical 
system that includes either an exponential or a UDM po- 
tential is autonomous and it admits extra Noether sym- 
metries (see previous section). Our aim is to derive the 
predicted time dependence of the main cosmological func- 
tions namely a(t) and (p(t) [and thus of H(t) and uu(t)] 
in the scalar field cosmology. 

4.1. Analytical solutions of the hyperbolic potential 

Inserting eq. (j5"Tj) into eqs. (|55|) . (|4"5)) the Lagrangian and 
the Hamiltonian respectively become 

L = \ (f - i 2 ) - ~ {^ix 2 - u 2 y 2 ) (58) 



After some algebra, we find that 

. ( I2E + UJ! . , n I2E + lo 1 ~\ 
6»i = In I - J sinh 9 2 + J smli 9 2 + 1 1 

We would like to remind the reader that the UDM dark 
energy model is recovered for uj\ = 2lo 2 . Obviously, we 
can easily prove that the concordance A-cosmology is a 
particular solution 4 of the current hyperbolic potential 
with e = 1, w\ = u>2 = wa and (6±, 9 2 ) = (0, 0). 

Finally it is interesting to mention that the UDM cos- 
mological model has been tested against the latest cos- 
mological data (SNIa and BAO) in Basilakos & Lukes 
72]. In this paper the authors discussed the evolution 
of matter perturbations as well as the spherical collapse 
model. They also compared the UDM scenario with the 
traditional A-cosmology and they found that the UDM 
scalar field model provides an overall (global and local) 
dynamics which is in a fair agreement with that of the A- 
cosmology although there are some differences especially 
at high redshifts. 



2 i^ix 2 



w 2 y 



(59) 



Technically speaking, in the new coordinate system our 
dynamical problem is described by two independent hy- 
perbolic oscillators and thus the system is fully inte- 
grable. In particular, utilizing the Euler-Lagrange equa- 
tions in the new coordinate system the corresponding 
equations of motion can be written as 

x — io\x = 
y - ^iV = 

where -v/cJi and ^Juj2 are the oscillators' "frequencies" 
with units of inverse of time. It is routine to perform the 
integration to find the analytical solutions: 

x (t) — sinh (y/uJit + 6*i ) 



y (*) = J 2E + UJl sinh {y^t + 2 ) 



0J 2 



a 3 (t) = -[x 2 (t)-y 2 (t)\ 
3 



Ul2 



[sinli (v^T< + 6i)- 



and 



>(i) = \ -f-arctanh (—\ 




arctanh 



2E + uii sinh (y/uJ^t + 6 2 ) 
ll>2 sinh (yfoJit + 9i) 



where 9i and 9 2 are the integration constants of the prob- 
lem. The constant 0\ is related to 9 2 because at the sin- 
gularity (t = 0), the scale factor has to be exactly zero. 



4.2. Analytical solutions of the exponential 
potential 

Now based on the exponential potential of eg. (157]). we 
extend the analytical solutions found by Russo [70( by ta- 
king into account the presence of non-relativistic matter 
(cold dark matter) , p m = lpx ^ 0. It is important to em- 
phasize that Russo [701 ] provided analytical solutions only 
in the context of quintessence DE (e = 1) with p m = 
while we solve analytically, for a first time (to our knowl- 
edge), the current dynamical problem by treating dark 
energy simultaneously either as quintessence or phantom 
with p m 7^ 0. Since for the current potential the (x, y) 
coordinate system does not lead to an analytical solution 
we are using the same methodology with that provided 
by Russo [70(| , we change variables from (a, </>) to (u,v) 
according to the transformations 



/^ + iln( a 3) 



(60) 
(61) 



Note that we have interchanged (it, v) with respect to 
those of [7(|. Inverting the above equations and using 
eq. ([57"l) we get 



a = e 3 




4 In the ACDM cosmology the scale factor is a(t) = 
ao sinn 2 / 3 (u>\t) where wa = 3 -Ho ao = (f^m/HA) 1 / 3 and 



On 



1 - fin 
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and 



4.2.1. Case K = 1 



V(u,v) = V e- d(u - v)/2 



(63) 



In the new variables (u, v) our Lagrangian ()33|) is writ- 
ten as follows: 



Ju+v) 



+ kV e- 2K ^ 



, K = - (64) 



where kVo = 8/3. The next step is to consider a change 
in the time coordinate as follows: 



dr 

~dt 



3kV 



-K{u-v) 



(65) 



which implies: 



= dv = dvdr = 3kVo c _ K(u _ v) 
dt dr dt V 4 



du du dr 
dt dr dt 



3kV ( 



0_ e -K(u-v) 



where u' = # and v' = Obviously using the latter 
transformations, eq. ([M|) and ea. (|55"|) the action given by 
eq. (l34l takes the form: 



J d 3 x dre^e' K ^ («V + 1) . (66) 



S = - 



Now varying the action we arrive at 

u" + {\-K)u' 2 - (1 + K) = (67) 
v" + (1 + K) v' 2 - (1 - K) = 0. (68) 

Obviously, in the latter equations the variables u, v de- 
couple. In these variables, the Hamiltonian of the system 
becomes: 



and (HU) be- 



In this case the system of eas. ([6T| . 
comes: 

u" -2 = 
v" + 2v' 2 = 



e 2v («V - 1) = E 

and the solution is: 

U (t) = T 2 + ClT 

« M = \ In (2c 3 r) 

where the constants are related by the constraint: 

E 



E = C3C1 c 3 



ci 



Replacing we find: 



u (r) = r 2 + cit 



We note that the solution depends on one arbitrary pa- 
rameter c\ 7^ 0. If we choose c\ = E then we have the 
solution: 

u (r) = r 2 + Et 
v(t) = ilu(2r) 



from these follows: 



a 3 ( r ) = V2^e r2+Er 



E = e (u+v) e -K(u-v )( U V-1) 



(69) 



where E 1 7^ (or p m 7^ 0). As we have already stated 
the above system of equations has been derived also by 
[70l] in the case of quintessence dark energy (e = +1), 
and it is solved only for E = (see appendix A). Here 
we prove that the same equations are valid also in the 
case of phantom dark energy in which the scalar field is 
imaginary, however the potential, the scale factor and the 
FRLW metric are real as they should. 

We conclude that the FRLW metric (using eqs. [65] and 
in the coordinates (t, x^) is: 



ds 2 



3kV 



"t^W 2 + a 2 { T )dx l dx l . (70) 



Bellow we provide analytical solutions for the two differ- 
ent cases. 



(T) = 4V3^ [2T2 + 2ST ~ ln(2T)] 



4.2.2. Case K 1 

In this case the solution of the system is: 
1 



2(K-l) 
1 



u(t) 

V{T) 2(1 + *) 
where C = f_T^ , 



■In 



Cc\ Sinn 2 (yJ\K 2 - T\r H 
C^c 2 Sinn 2 U/\K 2 - l|r 

sinw K > 1 
sinhw < K < 1 



(71) 
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and 9 1 being the phase constant. Without loosing the 
generality we can select d\ = 0. Also, the two constants 
of integration satisfy the condition: E = C3C1. Next we 
may choose c 3 = E 7^ and have the solution: 



u (r) = - 
v(t) 



1 



2(K-1 
1 



■In 



CE 2 Sinn 2 [y/\K 2 - l|r 



2(1 + if) 



In 



C" 1 Sinn 2 L/\ 



K 2 - II 



Finally, we transform this solution to the coordinates 
a(r) ,<p (t). Doing so we obtain 



a d ( T ) = C l^T\E\ — Sinn [^/\K 2 -1\t 



the scalar held potentials hierarchy, being the po- 
tentials that admit extra integrals of motion, and 
therefore appear to be promising candidates for de- 
scribing the physical properties of dark energy as 
well as extracting useful cosmological information. 
The existence of the new Noether integrals can be 
used to simplify the system of differential equations 
(equations of motion) as well as to determine the 
integrability of the system. 

• Based on the exponential and hyperbolic potentials 
we find that the main cosmological functions, such 
as the scale factor of the universe, the scalar field, 
the Hubble expansion rate and the metric of the 
FRLW space-time are provided analytically. 



<P(-r) = -~J — - 



C K*=i\E\ ^Sinn (\/\K 2 - l|r 



-ia Li- 



It is interesting to mention that in the case of < K < 
1 at late enough times we obtain that t <~ ln< and thus 
the scale factor evolves as a{t) cx t 2 / 3 ^ 1 ~ R2 . The current 
solution of the scale factor can provide, a recent cosmic 
acceleration (a(t) > 0) for K E (^, 1). 



In a future work we plan to apply the same approach 
30 to f(R) cosmological models. 

Acknowledgments. We would like to thank M. Plio- 
nis and L. Perivolaropoulos for useful comments and sug- 
gestions. 



Appendix A: Matter density versus Hamiltonian 



5. CONCLUSIONS 

In this paper we propose to use a theoretical model- 
independent criterion, based on first integrals of motion, 
usually named Noether symmetries in order to discrim- 
inate the dark energy (quintessence or phantom) mod- 
els within the context of scalar field FLRW cosmology. 
This is possible via the geometrical symmetries of the 
space-time in which both gravity and dark energy live. 
In particular, following the general methodology of [49| 
(see also the references therein), the Noether symmetries 
are computed for 9 distinct accelerating cosmological sce- 
narios that contain a homogeneous scalar field associated 
with different types of potentials. Note that the free pa- 
rameters of the dark energy models studied here can be 
constrained by using the current cosmological data. In 
particular, one has to perform a joint likelihood analy- 
sis utilizing for example the SNIa data Q, the shift pa- 
rameter of the Cosmic Microwave Background (CMB) [6] 
and the observed Baryonic Acoustic Oscillations (BAOs; 
[73j). Such an analysis is in progress and will be pub- 
lished elsewhere. 

The main results of the current paper can be summa- 
rized in the following statements (see sections 3.2 and 
4): 

• We verify that all the scalar field potentials, studied 
here, admit the trivial first integral, namely energy 
conservation as they should. 

• We find that the exponential and the unified dark 
matter potentials occupy an eminent position in 



We remind the reader that if the matter pressure is 
constant then the dynamical system described by the 
general Lagrangian of eg. (f3"3"|) is autonomous. Therefore, 

one can easily prove that p m = — P rn . Indeed, utiliz- 
ing eq.© and eq. (|35[) we have after some simple algebra 
that: 



K 3 k 



2 v 



V{<t>) 



ah 2 + K 3 a = \a 3 ^ 2 + a 3 ^p m + ^a 3 V (<jy) 
6 3 3 



6 3 



K 3 a 



'-a 3 [P„ 



or 



and thus: 



E = -ka 3 [P n 



ka 3 



p 



which satisfies eq. tfTTj) . Note that the inequality p m > 
points that E < 0. The case of non-relativistic matter 
P m = implies p m (a) = j^. 
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Appendix B: Solutions of Harmonic Oscillator 
coupling with a free particle 

Here we consider either oji ^ 0, u; 2 = or ui\ = 0, 
W2 ^ 0. For the latter case one may checks [46j]. In both 
cases the system is equivalent to a pair of two dynamical 
systems a simple harmonic oscillator and a free particle. 
The Lie point and the Noether symmetries of the current 
system can be found in [51| . As an example for lo\ ^ 0, 
cj 2 = equations (|58)l and (|59l) become 



UJ\X 



E = {u'v' - i) e ("+«) e -iC(«-«) 



(CI) 



Due to the fact that E — eq. (|C1[) takes the form uV = 
1. Thus we consider the following cases: 
Case 1: For K — 1 the solution is: 



u (f ) = r 2 



«(t) = -ln(2r) 



or 



and 



x — loix — y = . 

The solution of the above system is 

x (t) = sinh (^/uJit + 6i) 
y (*) = y/2E + to 1 t 



or 



a 3 (0 = -[x 2 (t)-y 2 (t)] 



- [sinh 2 (V^Ti + 6> x ) - (2E + cj x ) i 2 ] 
V2-E + uTi \ 




arctanh L 

3fee \ sinh ( ynUit + 0\ 



*(*) = 



Note that due to a(0) = we have d\ = 0. 



a 3 (t) = V2re r2 



■W4^[2r 2 - 1 n(2 T) ] 



Case 2: For if 7^ 1 the solution is: 
1 



u(t) 
v(t) 



(K - 1 
1 



hi 



(Vl^ 2 -i|r + 0i) 



or (for 6»i = 0) 

« 3 (r) = 



In 



Coss sJ\K 2 - ljr + 0i 



Coss (VAT 2 - It) 



SinnU/l^-llr 



Appendix C: Exponential Potential versus empty 
space p m = (or _B = 0) 

In this appendix we would like to give the reader the 
opportunity to appreciate the fact that our solutions pro- 
vided in section 4 can be viewed as an extension of those 
found by Russo [70| for the quintessence (e = 1) dark 
energy with p m = (or E = 0). For either quintessence 
or phantom dark energy, the system of equations (1671) . 
and (|oT))) which we have to solve is: 

u" + {1 - K) u' 2 -(1 + K) = 
v" + {l + K)v' 2 -(1-K) = 



0(r) 



1 



hi 



a 3 (r)Sinn Vl^ 2 -l|' 



2 V 3fce 

where the quantity Sinn is given by feal7ip and 



Cossw 



cosw K > 1 
coshw < K < 1. 



We point that for e = 1 the current solutions coincide 
(modulus some constants) those of [70| , In the case of 
phantom dark energy (e = —1) the scalar field is imagi- 
nary, however the potential, the scale factor and the met- 
ric of the space-time are real as they should. 
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